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Abstract. In a companion paper, we investigated the question of the spatial origin of the cosmic rays detected 
in the Solar neighborhood, in the case of standard sources located in the Galactic disk. There are some reasons 
to believe that there may also be a large number of sources located in the halo, for example if the Galactic dark 
matter is made of supersymmetric particles or if Primordial Black Holes are present. These exotic sources could 
enhance the p, d or positrons above the standard background, indicating the existence of new physics. The spatial 
distribution of these hypothetical sources, though an important ingredient to evaluate these exotic signals, is 
poorly known. The aim of this paper is to point out that this discussion should not be disconnected from that 
of the propagation properties in the Galaxy. More precisely, we determine the regions of the halo from which 
a significant fraction / of cosmic rays antiprotons and antideuterons detected in the Solar neighborhood were 
emitted (we refer t o these regions as /- volumes) , for different sets of propagation parameters consistent with B/C 
data, as derived in lMaurin et alj l|20o3). It is found that some of them lead to rather small /-volumes, indicating 
that the exotic cosmic rays could have a local origin (in particular for a small halo or a large Galactic convective 
wind), coming from the solar neighborhood or the Galactic center region. It is also found that the dark matter 
density enhancement (spike) due to the accretion around the central supermassive black hole gives a negligible 
contribution to the exotic charged particle signal on Earth. The case of electrons and positrons is also discussed. 
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1. Introduction 

A great amount of work has been done these last twenty 
years on the astrophysical signatures that could unravel 
new physics. In the eighties, there were great hopes that 
the antiproton signal, which showed an excess at an en- 
ergy of a few hundreds of MeV in the first balloon ex- 
periments, could be such a signature. However, this hope 
was swept away by the progress in measureme nts - see 
e.g. BESS (lOrito et all l2000t iMaeno et £illl200ll) or he at 
(jBeach et al.ll200l|) and caprice (|Boezio et all 1200 1|) at 
higher energy - and a better determinati on of the cos- 
mic r ay propagation parameters (see e.g. iMaurin et al.l 
1)2002(1 ). It was shown that the measured antiproton 
flux was indeed compatible w ith the sole secondary 
standard spallative production l|Bergstrom et al.l Il999bt 
iDonato et al.ll200l|l (see the first paper for a comprehen- 
sive historical discussion and a panel of references dealing 
with exotic antiproton production). 

iDonato et all l|2000() showed that the antideuteron 
(d) signal could lead to a clearer signature of SUSY. 
However, as discussed in many other studies on SUSY an- 
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determinacy in the dark matter distribution, as well as 
its possible dumpiness l|Bergstr6m et al.l Il999a|) . might 
severely change the conclusions. In contrast, the Hawking 
evaporation of Primordial Black Holes (PBH) could also 
yield a new source of cosmic rays (|Maki et al.lll99(j) . but 
the precise shape o f the dark matter in this case is not cru- 
cial (|Barrau et al.l[2002. 200"^ . Nevertheless, in the latter 
case, it was shown that even considering only the propaga- 
tion parameters giving a good fit to B/C data, the remain- 
ing degeneracy for example in the di ffusive halo height h as 
sizeable effects on the primary flux l|Barrau et al.ll2002() . 

Hence, at least two different phenomena can affect 
the conclusions reached in papers dealing with exotic flux 
calculations. The first one, related to the spatial distri- 
bution of SUSY sources , is usually thoroughly discussed 
(|Bergstrom et al.lll999b|) . but the second point - namely 
the influence of various propagation parameters - is gen- 
erally skipped, due to the simplicity of the propagation 
models used. The aim of the paper is not to compare the 



2 



Maurin & Taillet: Spatial Origin of Cosmic Rays in Diffusion Models 



predicted p, D fluxes with observations for different series 
of models, but rather to point out which characteristics of 
the models actually play a role, in order to give some phys- 
ical insights and milestones for studies specifically devoted 
to exotic flux evaluations. 

We apply the method described in iTaillet fc MaurinI 



(|2003(1 to determine the volumes from which a fraction / 
of cosmic rays reaching the Solar neighborhood were emit- 
ted, or cquivalcntly the volumes that contribute to the 
fraction / of the total flux detected in the Solar neighbor- 
hood. These volumes will be referred to as the /-volumes 
throughout the paper. 

We find that depending on the diffusion parame- 
ters (evaluated froin a systematic study of standard CR, 
iMaurin et al.l (|2002l) ) as well as on the source spatial dis- 
tribution, the spatial origin of cosmic rays may be quite lo- 
cal, the particles detected in the Solar neighborhood hav- 
ing mostly been created a few kpc away from the Solar 
neighborhood in some cases, or a few kpc away from the 
Galactic center in others. 

2. Evaluation of the /-volumes 

In a companion paper l)Taillet fc MaurinI 1200,'^ . we pre- 
sented a method to compute the region from which a cos- 
mic ray detected in the Solar neighborhood has a given 
probability of originating. This method was applied to 
standard sources located in the disk, and we now use 
it for (exotic) sources in the halo. A schematic view of 
our model is presented in Fig. El where the isothermal 
dark matter profile has been superimposed on the Galaxy 
to compare their typica l scales (the reader is referred to 
iTaillet fc MaurinI l(200.'^ for all the details concerning the 
model, such as the functional form of the galactic wind 
and the geometry of the box). The probability that a par- 
ticle detected in the Solar neighborhood was emitted from 
any finite volume V can be computed as 



r{v\ro} 



(1) 



where the source distribution w{rs) has been introduced 
and Nr^{ro) is the density in resulting from a point 
source located in Vs- In this paper, we are interested in 
determining /-volumes, i.e. volumes V(/) from which a 
given fraction / of cosmic rays detected in the Solar neigh- 
borhood were emitted. They are defined by 



v{vif)K} = f 



(2) 



Actually, even for a given value of /, there are many differ- 
ent volumes, delimited by different closed surfaces, fulfill- 
ing this condition. We focus on the smallest of these vol- 
umes, precisely delimited by an isodensity surface. Monte 
Carlo integration is then particularly well adapted to eval- 
uate the integrals in Eq. ^ In a typical run, ^ 10^ points 
are required to reach a < 0.5% convergence and the inte- 
gral is performed inside all isodensity surfaces at once, so 
that the /-volume defined by Eq. Q are simple to recover. 



2.1. Influence of the propagation parameters 

The quantity A^isl'^o) appearing in Eq. ^ is evaluated 
by solving the diffusion equation with a point-like source, 
in the geometry depicted in Fig. O Propagation is af- 
fected, at different levels, by three effects: escape, galactic 
wind and spallations. First, escape happens when a par- 
ticle reaches one of the boundaries of the diffusive vol- 
ume. As discussed in the companion paper, this limits the 
range from which cosmic rays can travel to the Solar neigh- 
borhood. It was also shown that the side boundary plays 
only a minor role, and one can assume that the box has 
an infinite radial extension. Second, a convective wind Vc 
directed out from the Galactic plane blows the charged 
nuclei away, so that it is more difficult to reach the plane 
from high z sources. Finally, spallations may happen when 
a nucleus crossing the thin disk interacts with the inter- 
stellar matter. The nuclei are then destructed at a rate 
Find = 2/inisM-^'-CTinci- A particlc emitted from a remote 
source is more affected by spallations as it is likely to have 
crossed the disk many times before reaching the Solar 
neighborhood. In the companion paper, this effect was 
shown to be important for heavy species created in the 
disk. Here, we focus on very light species, having smaller 
cross-sections, which are mostly created in the halo. They 
are affected by the wind in the whole halo, i.e. from the 
moment of their creation, whereas they are only affected 
by spallations when they cross the disk, which is less likely 
for halo sources than for disk sources. As a result, spal- 
lations play only a minor role in the present study (this 
effect is nevertheless included in our treatment). 

When these three effects are taken into account, the 
density in O due to a Dirac source <5(r — r^) can be com- 
puted. Because of the cylindrical symmetry present for 
an infinite disk, it is equivalent to consider a source term 
5{z — Zs)S(r — rs)/2TTrs, which leads to 
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give the order of magnitude of the typical distance over 
which the associated process affects propagation. In prac- 
tice, large values of R have been used in Eq. ||2Jl so that 
the hypothesis i? ^ oo is actually recovered. The effect 
of escape, wind and spallations arc compared in Fig. ^ 
which shows the shape of the isodensity surfaces for two 
values L = 5kpc (left panels) and L = 10 kpc (right 
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Fig. 1. Isodensity surfaces in the (z, r) plane for i = 5 kpc 
and L = 10 kpc (side boundary i? = 20 kpc). Inner con- 
tours correspond to d'P{rs, Zs\0) / d?r g = 0.01 kpc~'^ and 
the contours are spaced by a factor 1/4. From top to bot- 
tom, Tap = 3.17 kpc (no wind), no wind and no spalla- 
tions, Tw = 2.935 kpc (no spallations). These numbers 
correspond; respectively, for a reasonable choice of K{E) 
at 1 GeV, to cr^p ~ 100 mb and K ~ 20 kpc Myr^i. The 
additional thick line in each panel delimitates contours 
P(V(/)|0) =99%. 



panels), and for typical values of and rgp. The value 
Vsp = 3.17 kpc has been retained because it corresponds 
to the antidcutcron destruction cross section for a typical 
value of the diffusion coefficient K = 0.03 kpc^ Myr^^. In 
the upper panels, one can see the shrinking of the contours 
in the vicinity of the disk, due to the effect of spallations. 
The effect of the wind is rather to flatten the contours, as 
can be seen in the lower panels. The probability density 
also decreases more rapidly when convection or spallations 
are included than when diffusion alone is considered. As a 
result, the 99%-volumes are reduced, as indicated by the 
thick lines. It is thus of importance to use realistic values 
for K{E), L and Vc in order to give confident /-volumes 
for real situations (cTinoi is not a free parameter, it solely 
depends on the species we consider). 

To summarize the previous results about the ori- 
gin of exotic primaries in diffusion/convection/spallation 
models: i) the pure diffusive regime provides an upper 
limit that is strongly dependent on the halo size; ii) the 
Galactic wind lessens the /-volumes: either propagation 
is convection-dominated - in this case, the origin depends 
only on the value of L and r^, i.e. Vc and K{E) - or it 
is escape-dominated and the geometrical upper limit (sole 
dependence on L, not K) is recovered; iii) spallations also 
systematically lessen the /-volumes: the heavier the nu- 
cleus, the larger its destruction rate, the closer it comes 
from. However, as a particle created in the halo is less 
likely to cross the disk, this effect is negligible compared 



to the wind for > rgp. We show below that all these 
effects are more pronounced for annihilating SUSY than 
for evaporating PBH because the density profile /idm('', z) 
appears with a square. 

2.2. Dark matter distribution 

The dark matter distribution in our Galaxy is poorly 
known, and several dark matter profiles can be used. The 
first constraint is that the observed rotation curve of our 
Galaxy is almost flat beyond a few kpc from the center. 
For a spherical halo, it follows that the density decreases as 
l/r^ outside the central regions. In the inner regions, the 
situation is far from clear. Numerical simulations indicate 
that the central distribution of dark m atter is cuspy, wit h 
a dependence with 7 ~ 0.5 — 1.5 ((Ghez et al.lll998l) . 
but this seems to be i n contradiction with observations 
(|Binnev fc Evan£ll200l|) . 

In the absence of a clear answer to this problem, we 
use several profiles for the Dark Matter distribution, with 
the generic form 

'°"^''^^=IV^^J Ug + (Vr^ + z^)" J 

where spherical symmetry has been assumed. Numerical 
simulations point toward singular profiles with 7 = 1.5. 
a = 1.5, e = 1 and Rc = 33.2 kpc Ijvloore et aLlll999l) or 
7= 1, a = 1, e = 2 and Rc = 27.7 kpc ( l)Navarro et alJ 
Il996() . hereafter NFW). We also considered an isothermal 
profile with 7 = 0,a = 2,e = l and Rc — 3 kpc (the mod- 
ified isothermal profile would give very similar results). 

As already said, exotic SUSY particles (resp. PBH) are 
supposed to fill (resp. follow) the dark matter halo profile 
hriMir, z). However, the nature of the cosmic ray creation 
process is different in these two cases, leading to very dif- 
ferent effective source terms, i.e. different weight w{r, z) 
in Eq. For evaporating Primordial Black Holes, the 
particle production is proportional to the density of the 
objects wpBaif, z) oc /idm(''7^)- In contrast, the produc- 
tion term for supersymmetric particles is proportional to 
the square of the density because two dark matter parti- 
cles must be present for annihilation to occur. In this case 
wsusvir, z) oc hBuir, z)^. 

They are displayed in Fig. H both for SUSY and PBH 
weight (see above). The Moore and NFW profiles are 
singular at the Galactic center, so that the source term 
is much stronger there. The probability that a cosmic ray 
detected in the Solar neighborhood was emitted from this 
region is enhanced for these profiles. A crude estimate of 
this effect is obtained by a mere count of the effective 
(PBH or SUSY) source numbers in this critical region. For 
example, in the range [0 — 2] kpc, a Moore profile leads 
to an enhancement x2.7 for PBHs and x90 for SUSY 
annihilations, compared to the isothermal case. Stretching 
this interval decreases the enhancement factor, and for [0— 
4] kpc, it is respectively xl.5 and x25, and finally for [0 — 
8] kpc, the numbers are xl.l and x20. The enhancement 
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Fig. 2. Effective source weiglit (PBH or SUSY) for several 
profiles (see text). 
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Fig. 3. Schematic view of our Galaxy: diffusive and con- 
vective propagation plus spallations in the thin disk. 
Effective primary exotic sources follow either the dark 
matter profile or its square (isothermal profile is depicted). 

is far smaller for PBH than for SUSY particles. Notice 
that the upper limits on the P BH density derived from 
antiproton flux measurements in lBarrau et al ] ll2002l) were 
of the same order of magnitude for an isothermal halo and 
for cuspy halos. This result is definitively not transposable 
to the SUSY case. 

This is not the final word. The center of our Galaxy 
contains a supermassive black hole (SBH) of a few 10^ 
Mq. During its formation, it probably accreted the sur- 
rounding dark matter, leading to a local enhancement of 
the density. Gondolo & Silk (1999) (hereafter GS) found 



that if the SBH grows adiabatically in the center of the 
Galaxy, the cuspy profile (yo(r) oc r^^ with < 7 < 2) 
becomes spiky and p{r) oc with 2.25 < A < 2.5 in a 
region of a few parsecs around the black hole. The pres- 
ence of the spike would have dramatic consequences for 
several predictions of the signal from a nnihilating dark 
matt er particles, e.g. 7 and neutrinos (iGondolo fc Silk 



199S ) or synchrotron em ission of e+e pairs l|(5ondok 



2nOC ; iBertone et al.l l200lh . The signal coming from the 
direction of the Galactic Center is obtained by integrat- 
ing along the line of sight, and the contribution of the 
central region is very different with or without a spike. 
In the case of the isothermal profile, the central region 
(around the SBH) contributes at the level of ~ 10^*^ 
whereas this contribution is greate r than ~ 10^ for a 
Moore profile l|Gondolo fc Silkl|l999l) . However, these re- 
sults are expected to be overoptimistic, and it is doubtful 
that such a spike exists in our Galaxy, as indicated by 
a more careful dy namical modelling of the SBH growth 
(|Ullio et all 1200 ij) . These authors review several effects 
(adiabatic growth versus instantaneous growth, models 
with off-centered black holes) and recover some results 
that were known before the Gondolo fc Silk paper: only 
the peculiar case in which the SBH forms adiabatically at 
the exact center of the dark matter profile can lead to an 
enhan cement such as descr ibed in GS. Finally, in a recent 
studv. lMerritt etld] l)2002(l have observed that, taking into 
account the quite large probability that the Milky Way ex- 
perienced a major merger in its history, the ensuing dark 
matter profile and resulting annihilation fluxes could be 
several order of magnitudes smaller than obtained with 
dark matter profile not disturbed by a SBH. 

The points discussed above are mostly relevant for par- 
ticles travelling in straight lines. For charged particles, due 
to the diffusive nature of propagation, the probability to 
come from a sphere S of radius r = 10 pc around on 
the Galactic center (~ 8 kpc away) is Jg{d'P/cfirs)(Prs, 
which is < 10~^° {dV/cPrg is given for example in Fig.^. 
Due to the very narrow scale where the SBH may affect 
the distribution , even enhancement such as obtained in 
IGondolo fc Siii] l|l999ft - and which is not very realistic - 
cannot yield a significant contribution for charged parti- 
cles. Eventually, the dark matter profile remains of impor- 
tance (isothermal or cuspy). In the following, most results 
will be presented for the isothermal case, the influence of 
the cusp being discussed at the end. 



2.3. f -volumes for SUSY and PBH weights and 
different values of L 

We now have all the elements to compute the /-volumes, 
inserting the source distributions described above in 
Eq. (Q. The function entering the integral does not 
possess cylindrical symmetry, so that the full three- 
dimensional integral must be computed. We first neglect 
spallations and galactic wind to consider only the effect of 
L. This parameter is expected to play an important role. 



Maurin & Taillet: Spatial Origin of Cosmic Rays in Diffusion Models 



5 



as the charged particles created outside of the magnetic 
halo of o ur Galaxy do not penetra te inside it and are not 
detected ()Barrau et al.l2002ll2003|) . Fig.Hshows the 99%- 




Fig. 4. Contours 7'(Vsusy, pbhI-Rq) = 99% origin for 
L = 2, L = 5 and L = 10 kpc. /-volumes have been re- 
spectively evaluated with weight wpBii{r,z) cx /iDM(f, 
(solid lines) and wgusYir, z) cx h'^-^(r, z) (dashed lines) 
for R ^ L, but the result remain mostly unchanged us- 
ing R ~ 2Q kpc (but in case of a halo size i = 10 kpc 
that requires i? = 30 kpc). Upper panel: Vsusy(99%) and 
^^pbh(99%) in the Xg — j/s plane (zg = 0). Lower panel: 
same quantities but in the Xs — Zg plane {ys =0). In both 
panels, the dot marks the Galactic center, and denotes 
the Sun location (it is set to Rq = 8 kpc). 



volumes in the Galactic plane Zs = Q (upper panel), for 
the PBH and SUSY case. Their shape reflects the fact that 
the probability density has a maximum at = Rq , while 
the effective source distribution peaks at the Galactic cen- 
ter O. Because of the quadratic dependence on /idm('', z)^ 
Vsusy(99%) is smaller than Vpbh(99%). Three halo sizes 
are displayed (i = 2, 5 and 10 kpc): for larger halos, 
the surfaces are more deformed towards the Galactic cen- 
ter (the contribution of this region to the flux is larger), 
whereas they remain grossly unaffected in the anti-center 
direction. This effect is less pronounced in the case of a 
PBH-likc source distribution. The same contours arc also 
plotted for Us = in the lower panel. The shapes arc al- 
most maximally distorted towards rectangular contours. 
This is less and less pronounced, as either L is enhanced, 
or larger powers of /idm('', z) are chosen. 

The figures above show clearly that we are only sensi- 
tive to a well-defined region of the source distribution: first 
to the region which is embedded in the diffusive halo, and 
then, even within this region, to a sub-region between the 
Galactic center and the Solar neighborhood. These sub- 



regions represent a fraction of the total number of sources 
given by 



/v(99%) ^(^' ^) '^^^ 



/ w{r, z) (Pr 



where j w{r, z) (Pr is the total number of sources. It is 
also of interest to compare the number of sources located 
in the same sub-regions to the number of sources in the 
diffusive halo 



/v(99%) ^(^' ^) ^''^ 



where Vcyi is the volume of the diffusive halo. The cor- 
responding numbers are given in Tabled for various halo 
sizes. The fraction /'^^'(L) decreases with L, much faster 





PBH / SUSY 


PBH / SUSY 


L = 10 kpc 


~ 1. ~ 1. 


0.023 0.76 


L = 5 kpc 


0.70 0.85 


0.010 0.54 


L = 2 kpc 


0.31 0.60 


0.002 0.21 



Table 1. Fraction of the number of exotic primaries emit- 
ted in Upbh(99%) and Ususy(99%) for various L, com- 
pared to the total number of exotic primaries emitted ei- 
ther in the bounded geometry (halo size L and radial ex- 
tension i? = 20 kpc) - denoted /'^^'(L) -, or in the whole 
dark halo - denoted /'°*(L). 



for PBH than for SUSY. This can be understood as the 
number of contributors outside of the dark halo core ra- 
dius rapidly vanishes for SUSY particles (see Fig. 0). As 
regards the results for /'°*(L), we remark that this num- 
ber is particularly small for PBH, i.e. only a very small 
fraction of primordial halos distributed in the Galaxy con- 
tribute to the charged primary cosmic rays detected in the 
Solar neighborhood. 

Finally, it is also interesting to give the fraction of 
primaries that escape before reaching the Solar neighbor- 
hood. It is defined as 



/°^=(V) = l_/dotoct(-y^ ^ ^ 



J^wjr, z)iVeyi(r|j?Q) 
J^w{r,z)N^{r\RQ) d^r 



.(7) 



where Ncy\{r\RQ) and A^oo(^|^0) arc respectively related 
to the flux of particles detected at Rq, in the cylindri- 
cal geometry and in an unbounded space, from the same 
sources emitting from inside the volume V. Estimations 
for L = 10 kpc and L = 2 kpc are compiled in Tab.|21 The 
trends are conform to intuition. Forming greater fractions 
of the detected flux requires more distant sources, the lat- 
ter more easily escape through boundaries. For large diffu- 
sive halo L, the fraction that escape increases more quickly 
for PBH sources than for SUSY sources, whereas the con- 
verse is true for small halos. This is related to the fact 
that one has to compare the shape and typical extension 
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L = 10 kpc 


L = 2 kpc 


/^b'h : V(50-90-99%) 


40-55-64% 


45-75-88% 


fih%Y- V(50-90-99%) 


49-52-55% 


59-92-95% 



Table 2. Fraction of primaries f^^"^ emitted from the (50- 
90-99)%- volumes of the cylindrical geometry (see above) 
that escape through upper and lower boundary located 
at L = 10 kpc or L = 2 kpc (for PBH and SUSY effec- 
tive source distribution), before they can reach the Solar 
neighborhood. 



of the source distribution to the parameter L. The fraction 
of primaries which are emitted inside V but which never 
reach the solar neighborhood is actually greater than f^'^, 
as even in the case of diffusion in unbounded space, there 
are many trajectories which start in V and never reach 
the solar neighborhood (diffusion in three dimensions is a 
transient process). 



(L = 2 kpc, L = 6 kpc and L = 10 kpc). We empha- 
size that the values of all these parameters come from the 
study of standard sources of cosmic rays and do not de- 
pend on the exotic sources, which do not produce B nor 
C. We do not take reacceleration and energy losses into 
account in this work. These effects, though necessary to 
study the spectra of cosmic rays, are not so crucial here 
as they only amount to a redistribution of the cosmic rays 
at different energies. A particle detected at an energy of 
1 GcV/nuc was just created at a slightly different energy 
and its origin is not drastically different. 

The values of r^p and arc plotted in Fig. [S] for an- 
tiprotons and antideutcrons. The left panel shows that 
propagation is convection-dominated (r^ ^ 1) at low 
energy when large S values are considered and escape- 
dominated at all energies for small 6. Notice that although 
at a given d, the quantity r^jL is fairly independent of L, 
the origin is definitely not the same for L = 2 kpc as for 
L = 10 kpc. The right panel shows that spallation is not 



3. Realistic propagation parameters 

The previous section considered simplified diffusion situa- 
tions with a typical value K ^ 0.03 kpc^ Myr~^. Actually, 
K{E) is energy dependent, and more precisely, 

K{E) = A'o/37^^ 

Here, 6 is the diffusion slope and Kq the normalization 
of the diffusion coefficient. In a previous study (see Paper 
la, lb), we show that various combinations of parameters 
Kf), 6, diffusive halo height L and Galactic wind magni- 
tude Vc are equivalent, in the sense that they give a B/C 
spectrum that is consistent with the observations. In this 
section, we use these combinations to provide a realistic 
range of values for and rgp and to explore the conse- 
quences on the origin of exotic primary antiprotons and 
antideutcrons. The heavier antinuclei will not be consid- 
ered here, as it was shown by Chardonnet, Orloff & Salati 
(1997) that their formation is suppressed because of the 
low probability of coalescence of many antinucleons. 

3.1. Parameter range allowed 

To compute the parameters introduced in Eqs. 0, the 
spallation cross sections of antiprotons and antideutcrons 
are taken from the Particle Data Group^. In this work, 
we only consider spallation on pure hydrogen. It would 
be straightforward to take into account the spallations on 
the Helium component of the interstellar medium, but the 
effect is too small to be worth the complication. The four 
parameters Kq, S, L and Vc are taken from our compre- 
hen sive study of standa rd secondary to primary B/C ra- 
tio ijMaurin et al.ll200^ . Three values (two extremes and 
a medium value) have been retained for both the diffu- 
sion slopes {S = 0.35, 0.60 and 0.85) and the halo sizes 

^ http://pdg.lbl.gov/ 




1 10 10^ 1 10 10^ 



Ek (GeV/ruc) Ek (GeV/njc) 

Fig. 5. Left panel: evolution of as a function of kinetic 
energy per nucleus for primary p and d; from top to bot- 
tom, S = 0.35, S = 0.60 and 6 = 0.85. The parameter 
Tw/-/^ = Xw, as well as rgp/L = Xsp, are not very sensitive 
to the halo size L (for d, only L = 6 kpc is displayed) but 
and Tgp do. Right panel: Vsp/L as a function of Sfe/nuc 
for the same i5 values and for the halo size L = 6 kpc. The 
values of are different between p and d because they 
depend on rigidity (through K), i.e. on Z/A (it is 1 for p 
and 1/2 for d). For rgp, there is an additional strong de- 
pendence on the species because of the destruction cross 
sections. 

the dominant effect for the light nuclei considered here. 
Only for large diffusion slopes d and more particularly for 
antideutcrons this effect becomes sizeable and comparable 
to the diffusive escape. The comparison of the two panels 
shows that spallations are always less efficient than con- 
vective wind or boundary escape. Finally, whatever the 
value of (5, propagation is escape-dominated above a few 
tens of GeV/nuc and the origin of primary cosmic rays is 
solely dependent on the halo size. 
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3.2. Antiprotons and antideuterons 

We are now able to draw the /-volumes for the realistic 
propagation parameters being considered. We focus on the 
antideuteron signal as it seems to be the most promising 
species to look for in cosmic rays. An interstellar energy 
of 1 GeV/nuc is chosen; the nuclei that reach the detec- 
tor are solar modulated so that they are detected with 
a final energy of 400 — 800 MeV/nuc, where the signal 
is the more interesting. Table |3| summarizes the values of 
and at this energy for antideuterons. The situa- 





(kpc) 


5 = 0.35 


S = 0.60 


5 = 0.85 


L 


= 10 kpc 


Tw — 


oo 


8. 


2.9 








21. 


7.6 


3.5 


L 


= 6 kpc 


T'w — 


oo 


5.5 


2.1 






^sp — 


15.5 


5.5 


2.6 


L 


= 2 kpc 


/"w — 


oo 


2.1 


0.85 






l^sp — 


6. 


2.2 


1.05 



Table 3. rw and Vsp for three halo sizes L and three dif- 
fusion slopes S: these numbers are for 1 GeV/nuc (inter- 
stellar energy) antideuterons. 



tion is very different for small or large S. For small values 
(corresponding roughly to a Kolmogorov power spectrum 
S = 1/3), only spallations affect the propagation (Vc = 0, 

= oo) and this effect was shown to be weak; for large 
S - the value S = 0.85 is the one preferred in our B/C 
analysis (Paper lb) -, models are convection-dominated 
though Tgp and have about the same strength. 

Fig. displays 7'cyi(V(/)|C') = 99% for the values re- 
ported in Tab. 121 For S = 0.35 (external contours), the ge- 
ometrical (upper limit) contours are recovered. However, 
for larger 6 (internal contours), these contours shrink. All 
comments made in Fig. 01 as regards halo size, or SUSY 
and PBH behavior, remain valid. Actually, the diffusion 
coefficient slope 6, as L for the geometrical limit, is a key 
parameter to trace back the CR origin, because of the 
values of it implies, through Vc and Kq. 

It is also of interest to have a closer look at the first 
% that contribute to the fiux. As the /-volumes with 
/ < 50% correspond to isodensity contours that are quite 
insensitive to the boundaries (or to other effects) they 
present the axial symmetry around the Xs axis, so that 
a single cut through, e.g. the Xs — Zs plane, delivers all 
the information about their shape. Fig. [3 displays the /- 
volumes / =10-25-50-75% for L = 10 kpc. The difference 
observed in Fig. (71 between small (lower panels) and large 
S (upper panels) is readily explained: a la rge value of d 
also co rresponds to a large value of Kq (see iMaurin et all 
1)2003!) for details), so that a greater wind is needed in or- 
der to prevent from too many spallations occurring at low 
energy. The net result is that the wind blow the particles 
away, reducing the effective zone from where they come 
from. This is not the case for small S where the geomet- 
rical limit (pure diffusion) is almost reached. The conse- 



99% 

L=2 kpc 




1 1 1 i 




^ J 

■ -" SUSY ; 

— PBH ~ 

, , , 1 T 


20 

X. (kpc) 




1 1 1 1 J 




■ - Susy; 

— PBH ~\ 



L= 6 kpc 



-10 knc 



20 

X. (kpc) 



20 

K (kpc) 



.(kpc) 



Fig. 6. 99%-volumes for exotic primaries (no side bound- 
aries). Upper panels: cut in the Zg = kpc plane; lower 
panels: cut in the j/s = kpc plane. Left panels correspond 
to L = 2 kpc. middle panels to L ~ 6 kpc and right pan- 
els to L = 10 kpc. In each panel, we plot either the PBH 
case (solid lines) or the SUSY case (dotted fines). From 
external lines to internal lines correspond the values of the 
diffusion coefficient slope S — 0.35, S = 0.60, 5 ~ 0.85. 



sothermal profile 




. (kpc) 



.(kpc) 



Fig. 7. 7'cyi(V(/)|0) = 10% - 25% - 50% - 75% for L = 
10 kpc in the Zg = plane (except for the 75%-volumes, 
other /-volumes with / < 50% are not deformed by the 
boundaries so that they present symmetry around the Xs 
axis). Upper panels correspond to d = 0.85, and lower 
panels to (5 = 0.35. Both the PBH case (solid lines) and 
the SUSY case (dotted fines) are plotted. 



quences for indirect dark matter searches are important. 
In the case of an isothermal profile (left panels), the parti- 
cles created in the Galactic center have a small probability 
to reach the detector for large 5, whereas the converse is 
true for small 5. In the latter case, the predictions and 
the limits that can be put on a supersymmetric signal de- 
pend heavily on the central shape of the dark matter halo, 
which is precisely the part we know the least about. These 
contours for smaller diffusive halo sizes L have not been 
presented; they have a smaller extent, meaning that we are 
less sensitive to the distribution of dark matter far from 
the Solar neighborhood. As a result, the question of the 
dark matter density profile cusp is less crucial for small L. 
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Similar contours for the NWF profile are drawn in the 
right panels of Fig. |7| Combining information from the 
above surfaces to the relative enhancement of sources go- 
ing from the isothermal case to the cuspy case allows sev- 
eral complementary remarks: for small 6, about half the 
SUSY Cosmic rays come from the range [0 — 3] kpc. Thus, 
the ~ 50 enhancement factor on the production provided 
by the cusp translates directly into a factor x 50 in the de- 
tected flux. For PBH case, the origin is less localized and 
the enhancement factor is smaller, so that the net gain 
is more probably about 10 — 20%. For large 6, contours 
look like boxes encompassing both the Solar position and 
the Galactic center. In the SUSY case, the addition of a 
cusp strongly deforms the box towards the Galactic center, 
but it is not straightforward to estimate the enhancement 
without considering specific values for the diffusion pa- 
rameters. For PBH, the contours, and hence the flux, are 
not expected to be very sensitive to the parameters. This 
discussion is of less importance for small halo sizes. 

From the above discussion, it appears that the 
most important parameters are L and Vc/2K. The 
value S ~ 1/3 (Kolmogorov spectrum) corresponding 
to Vr = 0. has been preferred these last years (see 
e.g. IStrong fc Moskalenkol lll998tl). However our previ- 
ous studies l)Maurin et all l200li l2002l |2003|) show that 
large values of 6, and non-null values of Vc, are pre- 
fe rred. This trend is confi rmed the most recent results 
of IMoskalenko et al ] 1I2OO2I) who now tend to prefer S = 
0.42 — 0.52. To conclude, if the value of S happens to be 
large or more precisely if a strong Galactic wind is pre- 
ferred, the discussions about the dark matter profile, in- 
cluding about the existence of a spike, are not so crucial. 
If conversely 6 is small (no wind), the dark matter cusp 
as well as the exact location of the Solar system should be 
accurately known before exploring the SUSY parameter 
space. 

Finally, all the remarks made for antideuterons in the 
previous sections apply as well for antiprotons. According 
to Fig.|Sl for a given (5 at a given energy, the corresponding 

is about twice its antideuteron value. The resulting /- 
volumes are larger than those for antideuterons, but the 
conclusions remain the same. 



3.3. Electrons and positrons 



five l i fetime given bv l|Aharonian et alJl995l : lAtovan et alj 
ll995HBaltz fc Edsiolll999l) 



Tioss ~ 300 Myr x 



1 GeV 



E 



lAharonian et al ] l)l995(l and lAtovan et all l)l995|) showed 
that in that case, all boundaries have negligible effects 
on positrons and electrons above a few GeV, so that 
the characteristic distance travelled by these species is 
''loss ~ V Ktioss (random walk through the tangled mag- 
netic fields), or 



no 



1 kpc X 



1 GeV 



K 



E 



0.03 kpc^ Myr"^ 



The result is an exponential cutoff that depends on the 
energy, i.e. the probability density rea ds dVr^dld?r„ oc 
exp(— 7-s/rioss)/?'s (see also Sec.4.3 in iTaillet fc MaurirJ 
([2003)). In the case considered here of sources in the 
whole diffusive volume, the normalized probability den- 
sity is given by 



(8) 



It is quite different from the case of a s ource distribution 
locate d in the disk only (see Eq. (12) in ITaillet fc Maurii] 
(|2003|) '). The resulting /-volumes (spheres) are given by 



^rad(r < rii„|0) = 1-1 



?'lim 
^loss 



exp 



nil 



no 



(9) 



It means that sources that contribute to the fraction / = 
(50 — 90 — 99)% of the detected positrons emitted in the 
halo are located inside the sphere of radius riim « (1-7 — 
4.8— 6.6) xrioss- For the realistic values of K{E) used above 
(see also Paper lb), we compile in Tab.^jthe range covered 
by r\oss dX E ~ 7 GcV. Because of the very small scale 





5 = 0.35 


5 = 0.60 


5 = 0.85 


L = 10 kpc 


1. 


.65 


.48 


L = Q kpc 


.85 


.55 


.41 


L = 2 kpc 


.53 


.35 


.26 



Table 4. The quantity rioss (kpc) is given for three halo 
sizes L and three diffusion slopes 5 at the total energy 
E = l GeV. 



Exotic sources in the halo also emit electrons and 
positrons. Positrons are more promising to study super- 
symmetric signals as the background of standard positrons 
is much lower than electron's (e+/ (e+ + e~) < 0.1), being 
predominantl y seco ndary. Recently, the heat experiment 
( Coutu ct al1 ll999l) reported an excess at about 7 Ge V 
(see also the mass-91 experiment. iGrimani et ah! lj2002l) ). 

These particles are lighter than nuclei, so that they are 
subject to much stronger energy losses, due to synchrotron 
radiation and inverse Compton. This results in an effec- 



involved along with the exponential decrease, /-volumes 
for positrons are expected to be only slightly deformed by 
the dark matter distribution, except for small 5 and large 
L whose 99%- volumes extend up to ~ 7 kpc. 

It is possible now to make a few quantitativ e comments 
on the HEAT results and on the conclusion of lBaltz et alJ 
1I2OO2I) about this signal. They argued that, defining a 
boost factor related to the dumpiness of dark matter, one 
can accommodate with e"*" data without enhancing too 
much the antiproton signal. The point is that antiprotons 
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come from further than positrons, so that if a clump ex- 
ists close to us, its contribution of antiprotons is averaged 
over a larger zone than positrons. A comparison of Figs.El 
and numbers presented above gives a relative distance 



/ o 

- n- 



0.1 



for all reasonable stationary propagation models. 
However, considering large or small (5, the effect of the 
dumpiness factor is expected to be different in different 
propagation models. Hence, t he enhancement fa ctor for 
the antiproton signal used in iBaltz et al.l (|2002f) should 
also depend on the diffusion efficiency, i.e. combination 
of diffusion plus convection (that is not considered in 
the above reference). To summarize, the relation between 
SUSY positron and antiproton signals is not straightfor- 
ward, if the dark matter halo is clumpy. Thus it seems a 
hard task to combine constraints from these two different 
signals, unless they are obtained with the same analysis. 
Depending on 8 and L, their origin is more or less local, 
and the size of the clumps as well as the typical distance 
between the clumps may be of importance. 

4. Summary and conclusions 

This paper analyzes the spatial origin of exotic particles 
created from a dark matter profile. We presented the /- 
volumes inside which a given fraction of the cosmic rays 
detected in the Solar neighborhood were emitted. At high 
energy [E ^ \ GeV/nuc), the shape of the isodensity 
surfaces is set by the geometry of the diffusive halo, in 
particular on its height L, the influence of the side bound- 
ary at r = i? being small. We then showed that the /- 
volumes defined are smaller when spallations and convec- 
tion are taken into account, but in a very different way: 
for particles in the diffusive halo, the wind exponentially 
decreases the probability of reaching the Galactic plane, 
whereas spallations have about a null effect on the latter. 
The parameters L and 2Vc/-?4r indicate whether the prop- 
agation is convection or escape-dominated. In Table |S1 we 
summarize the parameters that act as a cut-off in vari- 
ous situations. Two source distribution for the isothermal 
dark matter profile were considered: production related to 
the density of the source (e.g. PBH evaporation), or pro- 
duction related to the square of the density of the sources 
(e.g. SUSY annihilation). The 99%-volumes are strongly 
stretched toward the Galactic center, corresponding to the 
maximum of the source distribution. This follows from the 
competition between the effective source which is maxi- 
mum at the Galactic center, and the probability density 
which steadily decreases from our position Rq to reach 
~ 10~^ — 10~^ kpc~^ for purely diffusive regime (or even 
less when convection is included) at the Galactic center. 
The fluxes in the Solar neighborhood are found to be far 
more sensitive to the dark matter profile in the SUSY case 
than in the PBH case. In both cases, the side boundary of 
the diffusive volume is observed to play a negligible role 
as long as i? > 20 — 30 kpc. 



As a last step, realistic propagation parameters were 
implemented, and the key parameters were found to be 
the halo size L and the diffusion slope 8 (actually 14/it^o)- 
For the species considered here (antiprotons and an- 
tideuterons), spallations always play a negligible role in 
the origin. It was found that this origin is far more local 
in case of large 8 and small L than in case of small 8 and 
large L. Moreover, the shape of the dark matter distri- 
bution near the Galactic center does not matter so much 
for the PBH case, whereas it may be crucial for SUSY 
annihilating particles. We emphasized that in any discus- 
sion of the annihilation signal in charged particles, the 
propagation parameter 8 or more precisely, the presence 
of a Galactic wind, should be considered, with the same 
importance of the parameter L or the choice of the dark 
matter profile. 

Two last points are worth noting. First, even though 
the work presented here does not allow a quantitative 
estimation of the effect of possible dumpiness of the 
dark matter halo (for SUSY annihilations), we observed 
that the comparison between the electron and antiproton 
SUSY signals should involve a careful inspection of the 
corresponding boost factors. Second, whereas the use of 
B/C-induced propagation parameters is justified for stan- 
dard antiprotons ( corres ponding /-surfaces can be seen in 
iTaillet fc Maiiriiil l|200,'^ ). there is no guarantee that these 
parameters are valid in the /-volumes depicted here. 
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Appendix A: Numerical evaluation of the point 
source solution in Bessel basis 

One needs to evaluate numerically point source solutions 
such as 



1 



-kKR 



ntponc 

E 

1=1 



Jo(C.r/i?) sinh[0(L-z)/E] 
C.Jf(G) cosh(C.L/i?) 



(A.l) 



In the above expression, (r, z) is the position of the 8 
source in polar coordinate and K is the radial extension of 
the Galaxy. 7V^^'(z) can be evaluated for each i and need 
to be summed till the ritionc-th order, which should for- 
mally tend to infinity. For evident reasons, ritronc is chosen 
to be the smallest possible with the constraint that the re- 
built series A^^'*''(r, z) has reached a good convergence. In 
the case of 8(r) point source, profiles are singular near 
the source and convergence of the series appear to be v ery 
slow. The ansatz depicted in iTaillet fc MaurinI l)2003j) is 
useless as soon as sources are outside the disk. We present 
below two alternatives to evaluate this sum. 
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Cut-off 


Escape-dominated Convection-dominated 
(Xw > 1) (Xw < 1) 


Losses-dominated 
(e-e+ > GeV) 


Halo 
Radial 


L L* ^ 3K/V = 3rw/2 
min(7l,3L) min(i?,3L*) 





Table 5. Summary of the typical cut-off in z and r directions beyond where a cosmic ray cannot originate. The sole 
parameters that determine these cut-offs are L (halo size) and/or Xw = fw/L ~ related to the convective wind Vc via 
Tw = 2K/Vc, or Tioss related to the effective life-time for positrons and electrons. 



A.l. Average value of the oscillating series with r 

In analogy with classical Fourier analysis, resummation of 
coefficients provide oscillating behavior around the true 
value. This can be understood if we recall that at the n- 
th order, the function added is oc Jq{x = C,nP)' P lying in 
[0—1], the argument of Jq takes values x = {Ci, C2, . . . Cn}, 
i.e. at the n-th order, the corrective function has n roots. 
Thus convergence will be more quickly reached if for a 
given order TT-cutoff, instead of evaluating Ng^^{r,z), one 
averages 

where r„_i, r„ and r„-|_i are ordered realizations of r. 
The sole condition is that the {r„}"=i' - "™'°ff belong to 
the grid r = {0, i?/(2ncutoff), 2i?/(2ncutoff), . . . , R}, 
i.e. 2ncutoff linear steps between and 1. To summarize, 
around the oscillating value, if the appropriate step is cho- 
sen, it ensures that the averaged two points are not both 
above or below the true value, and furthermore, that two 
opposite extrema of the oscillating function are averaged. 

A. 2. Step-like source: 9 function 

An alternative way is to consider solution from a step-like 
source, e.g. 0(a — r), in order to smooth the problematic 
behavior observed near the origin for the S source. With 
the suitable normalization in the source term, i.e. 

9{a-r) 
qe(r,z) = ^^(z) , 

and using the property J pJo{p)dp = pJi(p), it leads to a 
solution which is equivalent to the delta solution Ng^^{r, z) 
- Eq. ljA.l|) -, as long as the distance of the observer Xo 
from the source satisfies the relation To ^ a. The Besscl 
coefficients of 5 and 9 solutions are related through 

Arf(.„) = 2xMl^xiV/(z„). (A.2) 
(Qa/R) 

The acceleration of convergence can be understood as, in 
Eq. l|A.2p . the additional term behaves at least as 1/i (Ji 
is bounded and Q ~ Here a should be taken such as 
to verify a/R (with i? = 20 kpc for the Galaxy, one 
can safely take a ^ 10 pc). 

Thus, a 0-like source slightly underestimates the re- 
sult close to r = r^, but this zone corresponds to very 
small volumes that add a negligible contribution when one 



evaluates integrated probabilities. For practical purposes, 
both methods (average or 9 source) give the desired re- 
sults with about the same number of Bessel functions, i.e. 
"-cutoff ~ 100. 
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